NASA/TM-2013-217999 



An Approximate Solution and Master Curves 
for Buckling of Symmetrically Laminated 
Composite Cylinders 


Michael P. Nemeth 

Langley Research Center, Hampton, Virginia 


May 2013 


NASA STI Program ... in Profile 


Since its founding, NASA has been dedicated to the 
advancement of aeronautics and space science. The 
NASA scientific and technical information (STI) 
program plays a key part in helping NASA maintain 
this important role. 

The NASA STI program operates under the 
auspices of the Agency Chief Information Officer. 

It collects, organizes, provides for archiving, and 
disseminates NASA’s STI. The NASA STI 
program provides access to the NASA Aeronautics 
and Space Database and its public interface, the 
NASA Technical Report Server, thus providing one 
of the largest collections of aeronautical and space 
science STI in the world. Results are published in 
both non-NASA channels and by NASA in the 
NASA STI Report Series, which includes the 
following report types: 


• TECHNICAL PUBLICATION. Reports of 
completed research or a major significant phase 
of research that present the results of NASA 
Programs and include extensive data or 
theoretical analysis. Includes compilations of 
significant scientific and technical data and 
information deemed to be of continuing 
reference value. NASA counterpart of peer- 
reviewed formal professional papers, but 
having less stringent limitations on manuscript 
length and extent of graphic presentations. 

• TECHNICAL MEMORANDUM. Scientific 
and technical findings that are preliminary or of 
specialized interest, e.g., quick release reports, 
working papers, and bibliographies that contain 
minimal annotation. Does not contain extensive 
analysis. 

• CONTRACTOR REPORT. Scientific and 
technical findings by NASA-sponsored 
contractors and grantees. 


• CONFERENCE PUBLICATION. 

Collected papers from scientific and 
technical conferences, symposia, seminars, 
or other meetings sponsored or co- 
sponsored by NASA. 

• SPECIAL PUBLICATION. Scientific, 
technical, or historical information from 
NASA programs, projects, and missions, 
often concerned with subjects having 
substantial public interest. 

• TECHNICAL TRANSLATION. 
English-language translations of foreign 
scientific and technical material pertinent to 
NASA’s mission. 

Specialized services also include organizing 
and publishing research results, distributing 
specialized research announcements and feeds, 
providing information desk and personal search 
support, and enabling data exchange services. 

For more information about the NASA STI 
program, see the following: 

• Access the NASA STI program home page 
at http://www. sti. nasa. gov 

• E-mail your question to help@sti.nasa.gov 

• Fax your question to the NASA STI 
Information Desk at 443-757-5803 

• Phone the NASA STI Information Desk at 
443-757-5802 

• Write to: 

STI Information Desk 

NASA Center for AeroSpace Information 

7115 Standard Drive 

Hanover, MD 21076-1320 


NASA/TM-2013-217999 



An Approximate Solution and Master Curves 
for Buckling of Symmetrically Laminated 
Composite Cylinders 


Michael P. Nemeth 

Langley Research Center, Hampton, Virginia 


National Aeronautics and 
Space Administration 

Langley Research Center 
Hampton, Virginia 23681-2199 


May 2013 


Available from: 

NASA Center for AeroSpace Information 
7115 Standard Drive 
Hanover, MD 21076-1320 
443-757-5802 



Summary 


Nondimensional linear-bifurcation buckling equations for balanced, symmetrically laminated 
cylinders with negligible shell-wall anisotropies and subjected to uniform axial compression 
loads are presented. These equations are solved exactly for the practical case of simply supported 
ends. Nondimensional quantities are used to characterize the buckling behavior that consist of a 
stiffness-weighted length-to-radius parameter, a stiffness-weighted shell-thinness parameter, a 
shell-wall nonhomogeneity parameter, two orthotropy parameters, and a nondimensional 
buckling load. Ranges for the nondimensional parameters are established that encompass a wide 
range of laminated-wall constructions and numerous generic plots of nondimensional buckling 
load versus a stiffness-weighted length-to-radius ratio are presented for various combinations of 
the other parameters. These plots are expected to include many practical cases of interest to 
designers. Additionally, these plots show how the parameter values affect the distribution and size 
of the festoons forming each response curve and how they affect the attenuation of each response 
curve to the corresponding solution for an infinitely long cylinder. 

To aid in preliminary design studies, approximate formulas for the nondimensional buckling 
load are derived, and validated against the corresponding exact solution, that give the attenuated 
buckling response of an infinitely long cylinder in terms of the nondimensional parameters 
presented herein. A relatively small number of “master curves” are identified that give a 
nondimensional measure of the buckling load of an infinitely long cylinder as a function of the 
orthotropy and wall inhomogeneity parameters. These curves reduce greatly the complexity of the 
design-variable space as compared to representations that use dimensional quantities as design 
variables. As a result of their inherent simplicity, these “master curves” are anticipated to be 
useful in the ongoing development of buckling-design technology. 


Symbols 



nondimensional modal amplitudes defined by 
equations (13)-(15) 

An, A 12 , A 2 2, A 66 

orthotropic-cylinder membrane stiffnesses, lb/in. 

c x 

nondimensional coefficient defined by equations (44) 

Du, D 12 , D 22 , D 66 

orthotropic-cylinder bending stiffnesses, in-lb 

D m ( ),D C ( ), D b ( ) 

nondimensional linear differential operators defined by 
equations (3)-(5) 

E, G 

Young’s modulus of elasticity and shear modulus of an 
isotropic material, respectively, psi 

E 

reference value of Young’s modulus of elasticity (see 
equation (44a)), psi 
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lamina moduli (see Table 1), psi 
cylinder thickness (see figure 1), in. 


F F G 

J -'L’ J- T’ VJ LT 


h 

I 


L 

me {1,2, 3, ...} 

m cr 

ne{0, 1,2,3, ...} 

n cr 

N,N cr 

n x , n: 


P. 


Per 


R 


x,y 


x,x cr 


p 


nondimensional wall nonhomogeneity measure defined by 
equation (19) 

cylinder length (see figure 1), in. 

number of axial half waves in buckle pattern (see equation (13)) 
value of m at buckling 

number of circumferential waves in buckle pattern (see equation 
(13) and figure 2) 

value of n at buckling 

circumferential wave number defined by equation (20c) and 
value of N at buckling, respectively 

compressive axial prebuckling stress resultant and 
corresponding value at buckling, respectively 

nondimensional loading parameter defined by equation (12) 
and corresponding value at buckling, respectively 

cylinder radius (see figure 1), in. 

cylinder coordinate system (see figure 1), in. 

quantity defined by equation (28) and value at buckling, 
respectively 

nondimensional stiffness-weighted length-to-radius ratios 
defined by equations (8) and (6), respectively 

nondimensional orthotropy parameter defined by equation (9) 


8F nondimensional buckling-stress function (see equations (1) 

and (2)) 

6W nondimensional radial-displacement field at buckling (see 

equations (1) and (2)) 

ri = ^ nondimensional cylinder circumferential coordinate 

0 lamina fiber angle (see figure 1), degrees 
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X, X ct stiffness-weighted axial-buckle aspect ratio defined by 

equation (23) and value at buckling, respectively 

A, A cr nondimensional stiffness-weighted axial half-wave length of a 

buckle pattern defined by equation (18) and value at buckling, 
respectively 

p nondimensional orthotropy parameter defined by equations (7) 

v m laminate membrane Poisson's ratio defined by equation (10) 

v LT lamina major Poisson's ratio (see Table 1) 

| nondimensional cylinder axial coordinate 

p stiffness-weighted radius-to-thickness ratio defined by 

equation (11) 

a“ = value of the uniform axial stress at buckling, psi 


Introduction 

As evidenced by recent studies, efforts to develop shell-buckling design technology for 
laminated-composite cylinders are still underway. 1,2 Despite these and other efforts, the 
technology still has not reached the level of maturity needed to obtain the full weight-saving 
potential of shell structures, especially those made of laminated-composite materials. The major 
impediment is a lack of validated analysis and design methods. As a result, NASA SP 8007, 3 
developed in the late 1960s, is still the primary resource used by American industry for the design 
of buckling-resistant cylindrical shells. The procedure used in NASA SP 8007 has become a 
paradigm that accounts for the effects of initial geometric imperfections and edge support 
conditions by applying empirical “knockdown” factors to relatively simple classical linear- 
bifurcation buckling solutions. Although this class of solutions are very approximate, they 
generally provide insight into how variations in design parameters affect performance. Thus, it is 
useful to have relatively simple approximate analyses and formulas for predicting response trends 
that can be used as part of a building-block design process. 

Recent examples of design-technology development for buckling of laminated-composite 
cylinders have been presented in references 1 and 2. Specifically, baseline design formulas are 
presented for balanced, symmetrically laminated cylinders with walls that are effectively 
homogeneous and specially orthotropic. This approach was motivated by existing designs in 
which the allowable range of structural dimensions are such that the laminate is composed of a 
relatively large number of plies that are distributed uniformly through the thickness. Thus, it 
follows that modelling a cylinder wall as a homogeneous and specially orthotropic material is not 
a serious limitation in a preliminary design process, for a significant number of practical 
applications. Moreover, this approach is a logical first step in a building-block approach for more 
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general, laminate constructions with plies at the outer faces that are substantially different from 
those near the mid-plane. Most importantly, this approach yields significant simplification of the 
buckling equations and facilitates the development of simple design formulas in terms of useful 
nondimensional parameters, as evidenced by the results presented in references 1 and 2. Of 
particular importance are the “master curves” presented in these references that represent the 
behavior of a wide range of practical laminate constructions in a very concise manner. 
Specifically, the “master curves” presented in references 1 and 2 are plots of normalized buckling 
stress versus normalized buckling strain that are presented as a function of laminate composition. 
By including strain in these plots, both structural stiffness and load-capacity requirements can be 
considered in the design process. Moreover, curves such as these permit common behavior trends 
to be identified for a wide variety of laminated-cylinder configurations. The formulas and results 
presented in references 1 and 2 represent a significant first- approximation advancement in design 
technology for compression-loaded laminated-composite cylinders. 

The objective of the present study is to extend work presented in references 1 and 2 to include 
the cylinder-wall heterogeneity that cannot be neglected in a buckling analysis without 
introducing significant errors. Herein, cylinder-wall heterogeneity refers to laminate 
constructions in which there exists a significant through-the-thickness variation of the layer 
properties or orientations. Of particular interest herein are symmetric laminates with plies at the 
outer faces that are substantially different from those near the mid-plane. For convenience, this 
class of heterogeneous laminates are referred to herein as “sandwich-like” laminates. Use of the 
term “sandwich-like” herein is intended to connote the arrangement of plies and does not imply 
the existence of significant transverse shearing flexibility commonly associated with sandwich 
plates. To accomplish this objective, a set of nondimensional linear-bifurcation buckling 
equations and parameters are presented for balanced, symmetrically laminated cylinders 
subjected to uniform axial compression loads and in a membrane-stress prebuckling state. 
Anisotropy of the shell wall is presumed to be negligible. Then, a linear bifurcation buckling 
analysis is presented for the classical simply supported boundary conditions that are typically 
used in the preliminary design of cylinders. Next, approximate formulas for the buckling loads of 
infinitely long cylinders are derived. Results that establish the range of the nondimensional 
parameters appearing in the buckling analysis and simplified formulas are then presented. With 
the parameter ranges established, plots of nondimensional buckling load versus a stiffness- 
weighted length-to-radius ratio are presented for various combinations of the remaining 
parameters governing the behavior and for a very broad range of laminate constructions. Finally, 
relatively fewer “master curves” are presented for infinitely long cylinders that give a 
nondimensional measure of buckling as a function of the orthotropy parameters and the wall 
inhomogeneity measure. In addition, approximate formulas for the critical stresses are given that 
utilize the information in these curves. Thus, these “master curves” are anticipated to be useful in 
the ongoing development of buckling-design technology. 


Bifurcation-Buckling Equations 

To accomplish the objective of the present study, it is convenient to use the nondimensional 
parameters and equations presented in reference 4. When specialized to right circular cylindrical 
shells, such as that shown in figure 1, these equations correspond to Donnell-type, linear- 
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bifurcation buckling equations for symmetric laminates. In particular, the compatibility equation 
is given by 

D m (6F) = vT2 D C (6W) (1) 

and the transverse equilibrium equation is given by 

D b (6W) + vT2 D C (6F) = - (2) 


where D m ( ), D b ( ), and D c ( ) are nondimensional linear differential operators. The symbol 6W 
denotes the nondimensional radial-displacement buckling mode and 6F denotes the corresponding 
stress function. For symmetric laminates with negligible coupling between inplane stretching 
and shearing deformations, and between pure bending and twisting deformations, these operators 
are given by 


D m (6F) 


1 d 4 6F 

2 c _ 4 

a m d£ 


+ 2[i 


d 4 SF 

d^ 2 dr \ 2 


+ a„ 


d 6F 

dq 4 


(3) 


D C (6W) = p 


d 5W 


(4) 


D b (6W) 


i d 4 bw 

2 c 4 

a b dl§ 


+ 2p 


a 4 6w 

d§ 2 dq 2 


+ OCx 


d 6W 
dq 4 


(5) 


The nondimensional coordinates in these equations are defined as £ = A and q = ^ , where x is 


the axial coordinate, y is the circumferential arc-length coordinate, L is the cylinder length, R is 
the cylinder radius, 0 < % < 1 , and 0 < q < 2 jt. The nondimensional parameters appearing in these 
equations are given by 


a 


m 




( 6 ) 


F = 


A U A 22 ~ A 2 12 - 2A 12 A 

2A 66 \/A i i -A 2 2 


(7) 


a b 




( 8 ) 


n _ Pl2 + 2f^6 

vD||D 2 , 


(9) 
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V 


( 10 ) 


n/A n A 22 

1 

A u A 22 h . 

144D n D 22 J 1 

where A n , A 12 , A 22 , A 66 , D n , D 12 , D 22 , and D 66 are the stiffnesses of classical laminated-shell 
theory (see reference 5), and the corresponding stiffnesses A 16 , A 26 , D 16 , and D 26 are presumed 
negligible. The parameter defined by equation (11) is a stiffness-weighted shell-wall thinness 
metric. The quantity p is the nondimensional loading parameter defined by 


R / 

P =-yr~ 


_ N x R 2 (12) 

P 2 / 1 

Jt yD i ]b 22 

where N x is the magnitude of the uniform compressive axial stress resultant, that is presumed to 
exist prior to buckling, in classical linear bifurcation analyses for cylinders. 


Buckling Analysis 

For the purpose of preliminary design, cylinders with simply supported edges that restrict 
circumferential displacements are often used because they faciliate a simple analytical solution. 
These boundary conditions are usually referred to as the classical boundary conditions for simply 
supported shells and are denoted as “S2” boundary conditions in the recent book by Jones (see p. 
477 ) 6 nondimensional buckling modes for these boundary conditions are given by 


6W = ?4 sin(mjt§) cos(nr)) 


(13) 


where 0 < ^ < 1, 0 < r| < 2 jt, m G {1, 2, 3, ...} is the number of axial half-waves, n E {0, 1, 2, ...} 
is the number of full circumferential waves (see figure 2), and /f is the indeterminate modal 
amplitude. A typical buckling mode is shown in figure 3a that exhibits one full axial wave and 
four full circumferential waves. The corresponding contour plot of the radial displacement is 
shown in figure 3b. Buckling modes with n = 0 and n * 0 are referred to herein as 
axisymmetric and asymmetric modes, respectively. 


Substituting equation (13) into equation (1) indicates that the stress function can be represented 
by 


6F = S’ sin(m3t§) cos(nr|) 


(14) 


where S is related to the amplitude A. Specifically, substituting equations (13) and (14) into 
equation (1) yields 
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= - 


(15) 


2 S3 p m 2 jt 2 


4 4 

42 ~ 222 m Jt 

n a m + 2pn m jt + 




a n 


Substituting equation (15) into equation (14) and the result into equation (2), along with equation 
(13), and setting the coefficient of the trigonometric functions equal to zero yields the following 

equation for the eigenvalues p = p(m, n) 


4 2 2 2 

~ 2 n a b _ D 2 Jt m 
pjt = — 2 ~i + 2|3n + — 2 
m jt a„ 


+ 


12p 


4 2 

n a m 

2 2 

m jt 


+ 2jin + 


2 2 

2 m jt 


<X 


(16) 


For a given set of material properties and cylinder dimensions, the nondimensional parameters 
appearing in equation (15) are fixed values. The nondimensional buckling load, p cr , is then given 
by the smallest value of p that is found for m G {1, 2, 3, ...} and n G {0, 1, 2, The 
corresponding wave numbers are denoted by m cr and n cr . 


Approximate Buckling Formulas 

To obtain simple, approximate buckling formulas, similar to that given in reference 1, equation 
(16) is differentiated as if p = p(m, n) is a continuous function and m and n are real variables. 

Although not strictly correct, this approach gives useful results for certain ranges of the 
nondimensional parameters appearing in equation (16). To obtain a minimum set, 

d P = lE dm + lr dn = 0 (17) 


is enforced. In general, the partial derivatives appearing in equation (17) are very complicated 
functions of the variables m and n. To remedy this situation, it is convenient to introduce a 
nondimensional axial half-wave length A of a buckle pattern defined by 



(18) 


and a positive- valued wall nonhomogeneity measure given by the following ratio of extensional 
and bending stiffnesses 


( A i [1-^22 
V D n A 22 j 


(19) 
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For a homogeneous wall, 1 = 1. Because m is the number of axial half waves in a given 
bifurcation eigenvector, it follows that the quantity L/m is the axial half-wave length of the 
eigenvector. Substituting equations (18) and (19) into equations (6) and (8) yields 

(M=a < 20a > 

(jffif'Al (20b) 

As the cylinder length increases to infinity, A behaves like a continuous variable and 
minimization of p by differentiation with respect to A represents solutions for an infinitely long 

cylinder that buckles into a local shell mode and not into an overall column mode. Likewise, it is 
convenient to introduce 

N = n 2 (20c) 

By using equations (20), equation (16) becomes 


pjt 2 = N 2 AI + 2|3N + + 


12p 

2 1 
AN + 2pN + — 


( 21 ) 


Now, it is observed that application of the chain rule of differentiation to equation (17) reveals that 
minimization of p with respect to A and N is equivalent to minimization with respect to m 

and n, respectively. Enforcing = o and simplifying gives 


f 2 1 ) 

— + 2p, N + A N 2 

2 

2 

2 1 

IN 2 

- 12p 

N - — 

l IA j 

A 

A 


( 22 ) 


Likewise, enforcing = o gives 


(P + ian) 


-i- + 2pN + AN 2 
A 


- 12p 2 (u + AN) = 0 


(23) 


for values of N > 0, which yields 


-±- + 2|xN + AN 2 
A 
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12p 2 


U +AN 
P+IAN 


(24) 



Substituting equation (24) into equation (22) gives 

[i 2 a 2 n 2 - i](n + an) - i ( (3 + ian)[a 2 n 2 - i] = o (25a) 

for values of N > 0. For the case of N = 0, minimization of equation (21) with respect to nonzero 
values of A is given by equation (22), which reduces to 

+ 12p 2 = 0 (25b) 

IA 


Here, it is important to note that values of A and N that satisfy equations (24) and (25), denoted 
by A cr and N cr , respectively, yield the critical value of p , denoted by p cr , when substituted into 
equation (21). 

Solutions that minimize the right-hand side of equation (21) are obtained by examining the 
cases of N = 0 and N > 0 separately. For the case of N = 0, equation (25b) yields 


2^/31 p 


(26) 


Substituting equation (26) for A and N = 0 into equation (21) gives the result 


Per 


4^3 p 

jt 2 yr 


(27) 


where it is implied that N cr = 0. This solution corresponds to short-wave-length, axisymmetric 
buckle patterns. 

Solutions of equations (24) and (25a) for N > 0 are found in the present study by defining a 
positive-valued quantity X given by 

X = IAN (28) 


so that equation (25a) becomes 

(p - In)x 2 + (l - I 2 )x + l(p - ip) =0 


(29) 


equation (24) becomes 


f + 2 ^f 


12p ip + X 

I p + x 



(30) 
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and equation (21) becomes 


~ 2 
pJT = 


f 

X + 2(3 + 

V 



12p 


X I 

1 ^ X 


N 


(31) 


This transformation allows N to be expressed in terms of the nondimensional parameters and the 
new variable X. In particular, equation (30) is solved for N to get 


N = 2pX </3l ^ y + + X I 1 ' + 2pLXX + X2 ] 


(32) 


Equation (32) is now used to eliminate N from equation (31), which yields p a function of X; 
that is, equation (31) becomes 


P = 


2 \/3lp /Ip, + X 


X' + 2|3X +1 |3 + X 

+ • 


X + 2plX + I Ip + X 


(33) 


The critical values of p are obtained by substituting the positive- valued solutions to equation 
(29), denoted by X = X cr , into equation (33). In particular, the solutions to equation (29) are given 
by 


X = 


I - 1 


2(p - Ip) 


1 



4l(p - Ip)(lp - p) 

J. H” 9 

(x-0 


(34a) 


for P - ip * 0 and 1*1. The corresponding values of p„ and N cr are given by 


2\/3lp /IP + X C 


jt 


P + X C1 


x cr + 2pX cr +1 P + X c 

+ 


X cr + 2pix cr + I Ip + X c 


(34b) 


N cr = 2pX cr a / 3l 


Ip + x c 

P + X cr 


[x c 2 r + 2plX cr +I 2 ] ' 


(34c) 


Additionally, equation (28) yields 


A = 


X c , 

IN„ 


(34d) 
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These results correspond to short-wave-length, non-axisymmetric buckle patterns. It is important 
to note that as X cr + Ip -» 0, p cr -> oo. 


For, the case of N > 0 and 1 = 1 and either p - in * 0 or p - in = 0 , the solution to 
equation (29) is given by X cr = 1 and equation (34b) reduces to 


4V3p / 1 + P 
n 2 / y 1 + n 


(35a) 


which is identical to the corresponding result presented in reference 1. Additionally, equation 
(34c) reduces to 


N cr 


V3 p 

\J{ \ + P ) ( ' + ) 


(35b) 


and equation (34d) becomes 

= V(i + p)(i + n) 

V3 p 


(35c) 


For this case, p„ — > oo as p — * -1 and the corresponding axisymmetric buckling solution given 
by equations (26) and (27) govern the response. For cases with p - in = 0 and 1*1, equation 

(29) reduces to X = - p, which yields physically inadmissible solutions for positive values of (1, 

since X > 0 by definition. In addition, equations (34b)-(34d) become singular. Solutions for this 
case are obtained by re-examining equations (21)-(23). First, equation (23) yields 


II + 2pN + AN 2 ] - 12p 2 

A J 


(p + AN) = 0 


(36) 


which is satisfied by selecting 


1 2\ 2 12p 

— + 2pN + AN = — — 

A I 


(37) 


Substituting equation (37) into equation (22) gives 


I — — J— r — 0 

1 / A 


(38) 
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Enforcing 1^1 requires A' 1 to go to zero, which means that the buckle patterns have axial half 
waves that are on the order of the cylinder length. Next, substituting A" 1 = 0 back into equation 
(37) gives 


2^N + AN 2 = ± 2 



(39) 


For this equation, and equation (21), it is presumed that N -» 0 such that N 2 A and NA are finite 
valued. As a result, equation (21) with |3 = Ip yields the result 

4 -J3I p 

Pcr= — (40) 

JT 

The results producing equation (40) indicate that long- wave-length buckle patterns exist. In 
addition, it is important to note that the analysis presented by equation (22)-(35) is not valid for 
very long cylinders that exhibit buckle patterns with axial half waves that are on the order of the 
cylinder length. For the general case of long-wave-length buckle patterns, the parameter A is 
presumed to approach infinity and the parameter N is presumed to approach zero such that N 2 A 
and NA remain finite. Applying these presumptions to equation (21) gives 


~ 2 2 1ZU 

pit = N AI + -f- (41) 

N A 

Minimization of this equation with respect to N 2 A gives 

4v5i’ p 

Pc, = — (42) 

JT 

which is identical to equation (40), with no restrictions on the values of I, (3, or p. For this case, 
the buckle pattern is axisymmetric; that is, N cr = 0. Note that the vT appears in the numerator of 

equation (42) and in the denomenator of equation (27) for the short-wavelength axisymmetic 
buckle pattern. 


Results and Discussion 

In the present section, results that indicate the possible ranges of the nondimensional 
parameters and coefficients defined by equations (6)-(l 1) and (19) are presented first. Based on 
these ranges, generic buckling results obtained by minimization of equation (16) with respect to 
m and n are presented entirely in terms of the nondimensional parameters defined herein. These 
generic results are based on very wide parameter ranges and, as a result, encompass a wide range 
of geometries and laminate constructions and also illustrate a wide spectrum of behavioral trends. 
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Next, “master curves” for infinitely long cylinders are presented. In particular, two nearly 
identical sets of results are presented that were obtained by minimization of equation (16) with 
respect to m and n for extremely long cylinders and by using the approximate formulas derived 
herein. 


Parameter and Coefficient Ranges 

The nondimensional parameters and coefficients defined herein have been presented for a wide 
variety of balanced, symmetric laminates in reference 7 for the nine material systems given in 
Table 1. Two laminate families of particular interest that are discussed in this reference are the 
[(±45/0/90) p ] s and [(0/90/±45) p ] s quasi-isotropic laminates, and the [(+0/— 0) p ] s and [(-0/+0) p ] s 
angle-ply laminates. The results presented in reference 7 for the [(+0/— 0) p ] s and the [(-0/+0) p ] s 
angle-ply laminates indicate the following ranges: 

0.279 < [3 < 2.76 

-0.788 <; 5.43 


0.341 < (A 22 /A n ) < 2.93 


I = 1 


0.446 s 


A„A 22 h 
I44D m D, 


0.999 


for the nine material systems, regardless of the value of the stacking sequence index p. Likewise, 
the results for the [(±45/0/90) p ] s quasi-isotropic laminates indicate: 

1.03 |3 2.22 

0.926 <; (D 22 /D n ) 1/4 <; 1.00 


0.958 < Vl -vl 


A i t A 22 h 

144D„D 22 


1 

4 


< 1.09 


for the nine material systems and for p = 1,2, ..., 8. Moreover, for the [(0/90/±45) p ] s quasi- 
isotropic laminates the results indicate: 


0.266 s (3 0.968 
0.862 <(D 22 /D n ) 1/4 < 1.00 
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0.875 s \[\ -vl 


A„A 22 h J 

144D„D 22 


0.956 


for the nine material systems and p = 1,2, ..., 8. For both families of quasi-isotropic laminates, 
(A 22 /A n ) 1/4 = 1 and u =1 , regardless of the lamina material system and the number of plies. 
Thus, it follows that the heterogeneity parameter has the range 0.862 sis 1.00 for both families 
of quasi-isotropic laminates, regardless of the lamina material system and the number of plies. 

Two other laminate families considered herein are the [(±45/0 2 ) p ] s axially stiff laminates and 
the [(±45/90 2 ) p ] s circumferentially stiff laminates with integer values for the index p = 1, 2, ..., 8. 
Results obtained in the present study for the [(±45/0 2 ) p ] s laminates yielded: 

1.04 s (3 ^ 2.24 

0.531 < p< 0.987 

1.01 sis 1.72 


0.687 < (A 22 /A n ) 1/4 < 0.941 


0.901 < sf\ -vl 


A hA 22 h 
144D„D 22 


i 

4 


<0.988 


for the nine material systems and p = 1,2, ..., 8. For the [(±45/90 2 ) p ] s laminates, the same ranges 

i 

4 

were obtained. In contrast, the ranges 0.582 < I < 0.989, 
1.06 < (A 22 /A n ) 1/4 < 1.45 were obtained for the nine material systems and p = 1,2, ..., 8. 

Sandwich-like laminates with pairs of 0-deg and 90-deg plies lumped between pairs of angle 
plies were also considered in the present study. In particular, results were obtained in the present 
study for [(+0/-0) 2 /(O/9O) 4 ] s laminates that revealed: 

0.083 £0* 1.58 

1.22 < p < 16.89 

0.543 < I < 1.50 

0.845 <(A 22 /A n ) 1/4 < 1.18 


for (3, p, and \/T~ 


A hA 22 h 


144D,,D, 
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1.18 


0.988 < Vl - vL 


A,,A 22 h J 

144D„D 22 


for the nine material systems and for 0 <; 0 90 deg. Results were also obtained in the present 

study for 24-ply [90 p /0 q ] s laminates that indicate: 

0.065 < |3 < 0.697 

1.71 sps 17.85 

1.04 sis 3.42 

0.568 <;(A 22 /A n ) 1/4 < 1.76 


.829 s Vl-vL 


A ,, A 22 h J 

144D„D 22 


1.69 


for the nine material systems and for p = 1,2, ..., 11. Similarly, results for 24-ply [0 p /90 q ] s 

laminates that exhibit the same ranges for |3, p, (A 22 /A n ) 1/4 , and \fl - 

obtained. In contrast, the range of the heterogeneity parameter given by 0.292 sis 0.962 was 
obtained for the nine material systems and for p = 1, 2, ..., 1 1. 


A hA 22 h 


144D„D, 


were 


Generic Buckling Results 

The nondimensional buckling loads obtained by minimization of equation (16) with respect to 
m and n are shown in figures 4-75 as a function of the stiffness-weighted length-to-radius ratio 
a m defined by equation (6). Six curves are shown in each figure that correspond to values of the 
orthotropy parameter p = - 1,0, 1,5, 10 and 20. Figures 4-75 are also composed of three groups 
of related curves. In particular, figures 4-27 correspond to results for the stiffness-weighted 
thinness ratio p = 25, figures 28-51 correspond to results for p = 100, and figures 52-75 
correspond to results for p = 1000. Within each of these three groups, 24 curves are presented for 
combinations of the values of the orthotropy parameter (3 = 0, 1,2, and 3 and for values of the 
nonhomogeneity parameter I = 0.25, 0.5, 1, 1.5, 2, and 4. In computing these results, equations 

(6), (8), and (19) were used to obtain « b = vTa m . 

The results in figures 4-75 generally show a series of festoon curves in which the amplitudes of 
the festoons attenuate to a negligible magnitude as the stiffness-weighted length-to-radius ratio a m 
increases. Likewise, the nondimensional buckling loads forming a given curve attenuate to the 
constant value associated with the corresponding infinitely long cylinder that buckles into a local 
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shell mode, not an overall column mode. The results in figures 4-75 also generally show that the 
festoons making up the curves become less prominent as (3 decreases, as p increases, or as I 
increases for each group of curves corresponding to a given value of p. Additionally, the festoons 
become generally less pronounced for a group of curves as the value of p increases. 


“Master Curves” 

As evidenced by the parameter ranges presented herein, the numerous results appearing in 
figures 4-75 are applicable to a very broad range of symmetrically laminated cylinders. However, 
it is desirable to reduce these results into a compact form suitable for preliminary design that 
captures the essential trends and that provides adequate estimates of the buckling resistance. This 
reduction was performed in the present study by observing that each festoon curve appearing in 
figures 4-75 has a useful, practical lower bound that is given by the buckling load of the 
corresponding infinitely long cylinder. These lower-bound buckling loads were obtained in the 
present study for selected values of |3, p, p, and I by minimizing equation (16) with respect to 
m and n for relatively large values of a m . Moreover, the curves in figures 4-75 with the 
pronounced festoons required significantly larger values of a m than the curves that exhibit rapid 
attenuation of the buckling load with increasing values of a m . This process typically involves 
minimization over a very large set of values for m and n. As a result, the determination of a 
buckling load following this process involves a large number of computations. 

The approximate buckling loads for infinitely long cylinders that are given herein by the 
smallest positive value of p cr obtained from equations (27), (34b), and (42) require significantly 

fewer computations than the minimization of equation (16). To demonstrate the accuracy of these 
approximate buckling loads, several attempts were made during the present study to find a concise 
means of comparison applicable to the wide range of values for |3, p, p, and I appearing in 
figures 4-75. The comparisons presented in figures 76-79 were found to be adequate. An 
important step in arriving at the form of these figures was found when parametric studies revealed 

that plots of p cr versus p have the same shape for a fixed value of (3 and any value of p. In 
addition, it was found that identical sets of curves were produced for a given values of I and its 

lesf 

reciprocal. Moreover, equation (35a) revealed that - Pcr = \ for an isotropic material. These 

4</3 p 

2 r* 

facts led to the use of _ P cr as the ordinate in figures 76-79. 

4V3 p 

Four pairs of red and blue curves are shown in each of figures 76-79. Each pair corresponds to 
a given value of the orthotropy parameter (3. The blue curves correspond to buckling loads 
obtained by minimizing equation (16) with respect to m and n for a m = 100. The red curves 

correspond to the smallest positive value p„ obtained from equation (27) for short-wave-length, 

axisymmetric buckle patterns, from equation (34b) for short-wave-length, non-axisymmetric 
buckle patterns, and from equation (42) for long-wave-length, axisymmetric buckle patterns. All 
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of the curves shown in these figures were computed for a value of p = 100. In solving equation 
(34b), it was found that equation (34a) exhibited only one positive-valued root. 

Inspection of the results used to obtain figures 76-79 indicated differences in buckling loads for 
each pair of curves that are less than 3%. Similar results were obtained for values of p = 1000 
and 25 which exhibited differences in buckling loads for each corresponding pair of curves that 
are less than 3% and 8%, respectively. Overall, the differences in each pair of corresponding 
buckling loads obtained from the two solution methods are insignificant compared to the inherrent 
variability in the material properties. It is important to observe that the results in figures 76-79 are 
applicable to an enormous number of laminate constructions. Moreover, the compact form of 
figures 76-79 make the design space easier to navigate and interpret. 


The results presented in figures 76-79 show that the nondimensional buckling load n Eg 

4 \/3 p 

decreases as the orthotropy parameters (3 or u decrease. Moreover, the effect of decreasing u 
on the nondimensional buckling load is reduced as the nonhomogeneity parameter I becomes 
greater or less than a value equal to one. The further I departs from 1 = 1, the more the 
nondimensional buckling loads are reduced. 


Approximate Formulas for Critical Stresses 

The results presented in figures 76-79 are also useful for estimating the critical stress in a 
cylinder incipient to buckling. Specifically, approximate formulas for the critical stress of a 

geometrically perfect cylinder, a" are obtained by using a" = with equation (12) to get 


~ 2 
cr _ p cr TC 
° x — 2 

R h 




(43) 


An insightful, alternate form of equation (43) is obtained by using equation (1 1) to get 


where 


cr 




(44a) 


c = 


1 - v r 


144 A n A 22 D n D 2 
E h 


2 ~ 

rc p cr 
4V3 p 


(44b) 
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and the term in parentheses is the ordinate used in figures 76-79. The symbol E denotes a 
convenient reference modulus of elasticity and is included in equations (44) to provide some 
freedom in the way one may want to represent the critical stress for a given problem. For 
example, consider an isotropic material with a modulus E and Poisson’s ratio v. For this 
2 ^ 

material, - P° r = | and using E = E in equations (44) yields the well-known formula 
4 s/3 p 


Ox = 


V 3 ! 1 - y2 ) 


Eh 

~R 



( 45 ) 


that is often used in the design of isotropic cylinders. For a given class of laminated cylinders, 
picking the reference modulus to be the effective axial modulus of a specific laminate configura- 
tion in the class may yield a desired physical interpretation of equation (44a) that is useful in the 
design of these cylinders. 


Concluding Remarks 

A set of nondimensional linear-bifurcation buckling equations for balanced, symmetrically 
laminated cylinders subjected to uniform axial compression loads have been presented that are 
based on Donnell’s equations. These equations presume negligible shell- wall anisotropies and 
have been solved exactly for the practical case of simply supported ends. Nondimensional 
quantities have been used to characterize the buckling behavior that consist of a stiffness- 
weighted length-to-radius parameter, a stiffness-weighted shell-thinness parameter, a shell-wall 
nonhomogeneity parameter, two parameters that quantify membrane and bending orthotropies, 
and a nondimensional buckling load. Results that establish the range of the nondimensional 
parameters appearing in the buckling analysis have also been presented that encompass a wide 
range of laminated-wall constructions. Numerous generic plots of nondimensional buckling load 
versus a stiffness-weighted length-to-radius ratio are presented for various combinations of the 
remaining parameters governing the behavior. These plots span the established ranges of the 
nondimensional parameters and are expected to include many practical cases of interest to 
designers. Additionally, these plots show how the parameter values affect the distribution and size 
of the festoons forming each response curve and how they affect the attenuation of each response 
curve to the corresponding solution for an infinitely long cylinder. 

To aid in preliminary design studies, approximate formulas for the nondimensional buckling 
load have been derived. These formulas give the attenuated buckling response of an infinitely 
long cylinder in terms of the previously described nondimensional parameters and have been 
validated against the corresponding exact solution. An important finding of the present study is 
that a relatively small number of “master curves” have been identified that give a nondimensional 
measure of the buckling load of an infinitely long cylinder as a function of the orthotropy and wall 
inhomogeneity parameters. In particular, these curves reduce greatly the complexity of the 
design-variable space as compared to representations that use dimensional quantities as design 
variables. Approximate formulas for the critical stresses have also been given that utilize the 
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information in these curves. As a result of their inherent simplicity, these “master curves” are 
anticipated to be useful in the ongoing development of buckling-design technology. 
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Table 1. Lamina Properties 


Lamina 

property* 

Material Systems 

Boron- A1 

S-glass- 

epoxy 

Kevlar 

49-epoxy 

IM7/ 

5260 

AS4/ 

3502 

AS4/ 

3501-6 

Boron- 

epoxy 

IM7/ 

PETI-5 

P-100/ 

3502 

E l , Msi 

33 

7.5 

11.02 

22.1 

18.5 

20.01 

29.58 

20.35 

53.5 

E v Msi 

21 

1.7 

0.8 

1.457 

1.64 

1.30 

2.68 

1.16 

0.73 

V LT 

0.23 

0.25 

0.34 

0.258 

0.30 

0.30 

0.23 

0.29 

0.31 

G lp Msi 

7.0 

0.80 

0.33 

0.860 

0.87 

1.03 

0.81 

0.61 

0.76 


* The symbols L and T denote the longitudinal fiber and transverse matrix directions of a specially orthotropic lamina, respectively. 



Figure 1. Cylinder geometry, surface coordinate system, and lamina fiber orientation. 
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Figure 2. Circumferential waveforms used in equation (12). 
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(a) Three-dimensional rendering of a buckling mode 


23tR 


Axial half-wave length 



(b) Contour plot of the radial displacement 


Figure 3. Typical buckle pattern of a compression-loaded cylinder. 


22 


50 


Per ” 


. , cr 2 

Nx R 

n VOhD 2 


30 


20 


10 


ji = -1 


p = K V! -v 2 „ 


J = J ^11^22 

Id„A 2 2, 


144D„D 22 


p _ D12 + 2 D 6i 


^ _ A 11 A 22 A 12 2A 12 A 6 

2A 66 A , , A 22 


.0 




1 

/5 













^10 

O 

CM 

II 


= jjA 


Stiffness-weighted length-to-radius ratio, « m = r a 


Figure 4. Effects of parameters a,,, and p. on nondimensional buckling loads on 
cylinders with p = 0, 1 = 0.25, and p = 25. 
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Figure 5. Effects of parameters a m and \i on nondimensional buckling loads on 
cylinders with P = 1, 1 = 0.25, and p = 25. 
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Figure 6. Effects of parameters a m and \i on nondimensional buckling loads on 
cylinders with P = 2, 1 = 0.25, and p = 25. 



Stiffness-weighted length-to-radius ratio, a m = 


Figure 7. Effects of parameters a m and \i on nondimensional buckling loads on 
cylinders with P = 3, 1 = 0.25, and p = 25. 
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Stiffness-weighted length-to-radius ratio, « m = p(^ 22 ] 


Figure 8. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with P = 0, 1 = 0.5, and p = 25. 



Figure 9. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with P = 1, 1 = 0.5, and p = 25. 
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Figure 10. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with |3 = 2, 1 = 0.5, and p = 25. 
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Figure 11. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with P = 3, 1 = 0.5, and p = 25. 
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Figure 12. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with p = 0, 1 = 1, and p = 25. 
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Figure 13. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with (3 = 1, 1 = 1, and p = 25. 
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Stiffness-weighted length-to-radius ratio, = j^f^J 


Figure 14. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with P = 2, 1 = 1, and p = 25. 



Stiffness-weighted length-to-radius ratio, a m = j^f^J 


Figure 15. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with P = 3, 1 = 1, and p = 25. 
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Figure 16. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with P = 0, 1 = 1.5, and p = 25. 



Stiffness-weighted length-to-radius ratio, a m = j^f^J 


Figure 17. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with p = 1, 1 = 1.5, and p = 25. 
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Stiffness-weighted length-to-radius ratio, 


Figure 18. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with |3 = 2, 1 = 1.5, and p = 25. 
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Figure 19. Effects of parameters a m and g. on nondimensional buckling loads on 
cylinders with P = 3, 1 = 1.5, and p = 25. 
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Stiffness-weighted length-to-radius ratio, 




Figure 20. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with p = 0, 1 = 2, and p = 25. 



Figure 21. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with P = 1, 1 = 2, and p = 25. 
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Figure 22. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with P = 2, 1 = 2, and p = 25. 
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Figure 23. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with P = 3, 1 = 2, and p = 25. 
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Figure 24. Effects of parameters a,,, and u on nondimensional buckling loads on 
cylinders with p = 0, 1 = 4, and p = 25. 


Per “ 


20 


N! r R 2 15 


n VOhD 2 


10 


p 

AiiA 22 h 

i 

P _ D 12 + 2D 66 

p = ^Vi - v 2 m 

144D„D 22 

\/D„D 22 

T-( A u D 22 f 

A„ 

II — 

iA 22 - A 12 — 2A 12 A 66 

iDuAaj 

M" — 

2A 66 VAhA 22 


u = -1, C 

1, and 1 

















\ 




10 

O 

CM 

II 


= jjA 


Stiffness-weighted length-to-radius ratio, « m = ^ — 


Figure 25. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with P = 1, 1 = 4, and p = 25. 
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Figure 26. Effects of parameters a,,, and u on nondimensional buckling loads on 
cylinders with P = 2, 1 = 4, and p = 25. 
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Figure 27. Effects of parameters a m and p. on nondimensional buckling loads on 
cylinders with p = 3, 1 = 4, and p = 25. 
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Stiffness-weighted length-to-radius ratio, a m = 


Figure 28. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with P = 0, 1 = 0.25, and p = 100. 
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Figure 29. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with P = 1, 1 = 0.25, and p = 100. 
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Stiffness-weighted length-to-radius ratio, a m = 


Figure 30. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with P = 2, 1 = 0.25, and p = 100. 
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Figure 31. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with P = 3, 1 = 0.25, and p = 100. 
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Stiffness-weighted length-to-radius ratio, a m = 


Figure 32. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with p = 0, 1 = 0.5, and p = 100. 



Figure 33. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with (3 = 1, 1 = 0.5, and p = 100. 
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Figure 34. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with P = 2, 1 = 0.5, and p = 100. 



Figure 35. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with P = 3, 1 = 0.5, and p = 100. 
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Stiffness-weighted length-to-radius ratio, a m = 


Figure 36. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with P = 0, 1 = 1, and p = 100. 
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Figure 37. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with P = 1, 1 = 1, and p = 100. 
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Figure 38. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with p = 2, 1 = 1, and p = 100. 
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Figure 39. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with P = 3, 1 = 1, and p = 100. 
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Figure 40. Effects of parameters a,,, and u on nondimensional buckling loads on 
cylinders with p = 0, 1 = 1.5, and p = 100. 
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Figure 41. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with (3 = 1, 1 = 1.5, and p = 100. 
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Stiffness-weighted length-to-radius ratio, a m = 


Figure 42. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with P = 2, 1 = 1.5, and p = 100. 
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Figure 43. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with P = 3, 1 = 1.5, and p = 100. 
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Figure 44. Effects of parameters a,,, and pi on nondimensional buckling loads on 
cylinders with p = 0, 1 = 2, and p = 100. 
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Figure 45. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with P = 1, 1 = 2, and p = 100. 
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Figure 46. Effects of parameters a,,, and u on nondimensional buckling loads on 
cylinders with p = 2, 1 = 2, and p = 100. 


_ ,cr _ 2 

N v R 


Per 2 n n 

Jl V^11^22 


50 


40 


30 


20 


10 


£ 


p = -1, 0, and 1 


V. 


R 

AhA 22 h 

1 

o _ D 12 + 2D 66 

P = K Vl-v m 

144D„D 22 

Vd„d 22 

>-|CM 

CM 

CM 

o 

h 

A„ 

11“ 

CD 

CD 

< 

CM 

< 

CM 

1 

CM 

V 

1 

CM 

CM 

< 

V ^11^22 J 

M* — 

2 A 66 V A -1 -1 A 22 




10 


|i = 20 


= ijA 


Stiffness-weighted length-to-radius ratio, <* m = A 22 


Figure 47. Effects of parameters a m and \i on nondimensional buckling loads on 
cylinders with |3 = 3, 1 = 2, and p = 100. 
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Figure 48. Effects of parameters a,,, and u on nondimensional buckling loads on 
cylinders with (3 = 0, 1 = 4, and p = 100. 
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Figure 49. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with P = 1, 1 = 4, and p = 100. 
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Figure 50. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with P = 2, 1 = 4, and p = 100. 
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Figure 51. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with 0 = 3, 1 = 4, and p = 100. 
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Figure 52. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with P = 0, 1 = .25, and p = 1000. 
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Figure 53. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with P = 1, 1 = .25, and p = 1000. 
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Figure 54. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with P = 2, 1 = .25, and p = 1000. 



Figure 55. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with P = 3, 1 = .25, and p = 1000. 
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Figure 56. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with p = 0, 1 = 0.5, and p = 1000. 
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Figure 57. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with (3 = 1, 1 = 0.5, and p = 1000. 
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Figure 58. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with P = 2, 1 = 0.5, and p = 1000. 
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Figure 59. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with P = 3, 1 = 0.5, and p = 1000. 
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Figure 60. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with p = 0, 1 = 1, and p = 1000. 
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Figure 61. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with (3 = 1, 1 = 1, and p = 1000. 
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Figure 62. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with P = 2, 1 = 1, and p = 1000. 
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Figure 63. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with P = 3, 1 = 1, and p = 1000. 
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Figure 64. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with p = 0, 1 = 1.5, and p = 1000. 
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Figure 65. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with (3 = 1, 1 = 1.5, and p = 1000. 
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Figure 66. Effects of parameters a m and g. on nondimensional buckling loads on 
cylinders with P = 2, 1 = 1.5, and p = 1000. 
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Figure 67. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with P = 3, 1 = 1.5, and p = 1000. 
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Figure 68. Effects of parameters a m and g. on nondimensional buckling loads on 
cylinders with p = 0, 1 = 2, and p = 1000. 
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Figure 69. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with P = 1, 1 = 2, and p = 1000. 
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Figure 70. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with P = 2, 1 = 2, and p = 1000. 
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Figure 71. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with p = 3, 1 = 2, and p = 1000. 
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Figure 72. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with p = 0, 1 = 4, and p = 1000. 
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Figure 73. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with P = 1, 1 = 4, and p = 1000. 
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Figure 74. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with p = 2, 1 = 4, and p = 1000. 


Per = 


N 


cr 

x 


R 


2 


2 

jt 


\Dn D 22 



Stiffness-weighted length-to-radius ratio, a m = p(^ 22 | 


Figure 75. Effects of parameters a m and p on nondimensional buckling loads on 
cylinders with P = 3, 1 = 4, and p = 1000. 
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Figure 76. Nondimensional buckling loads of infinitely long cylinders with 1 = 1. 
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Figure 77. Nondimensional buckling loads of infinitely long cylinders with I = 0.5 and 1 = 2. 
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Figure 78. Nondimensional buckling loads of infinitely long cylinders with I = 0.67 and I = 1.5. 
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Figure 79. Nondimensional buckling loads of infinitely long cylinders with I = 0.25 and 1 = 4. 
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